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Abstract 

We study the tunneling of Bogoliubov excitations through a barrier in a Bose-Einstein conden- 
sate. We extend our previous work [Phys. Rev. A 78, 013628 (2008)] to the case when condensate 
densities are different between the left and right of the barrier potential. In the framework of 
the Bogoliubov mean-field theory, we calculate the transmission probability and phase shift, as 
well as the energy flux and quasiparticle current carried by Bogoliubov excitations. We find that 
Bogoliubov phonons twist the condensate phase due to a back-reaction effect, which induces the 
Josephson supercurrent. While the total current given by the sum of quasiparticle current and 
induced supercurrent is conserved, the quasiparticle current flowing through the barrier potential 
is shown to be remarkably enhanced in the low energy region. When the condensate densities are 
different between the left and right of the barrier, the excess quasiparticle current, as well as the 
induced supercurrent, remains finite far away from the barrier. We also consider the tunneling of 
excitations and atoms through the boundary between the normal and superfluid regions. We show 
that supercurrent can be generated inside the condensate by injecting free atoms from outside. 
On the other hand, atoms are emitted when the Bogoliubov phonons propagate toward the phase 
boundary from the superfluid region. 

PACS numbers: 03.75.Kk,03.75.Lm,67.85.De 
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I. INTRODUCTION 

In the pioneering work by Bogoliubov {l| , it was shown that the Bose-Einstein condensate 
(EEC) of weakly interacting bosons has a phonon-type excitation mode. It is now called 
the Bogoliubov mode, which is a Nambu-Goldstone mode associated with a spontaneous 
broken U(l) symmetry [2]. This collective mode dominates low-energy properties of BEG, 
so that it is an important key to understand physical properties of BEG [3]. In particular, 
the existence of Bogoliubov phonon is essential for the Bose-condensed_phase to acquire 

n r n 

superfluidity [J] . Since the realization of BEGs in ultracold atomic gases [5|, l6|l , the study of 
Bogoliubov mode has been one of the main issues in cold atom physics |3|, l7(]. Because of 
the high degree of controllability, the BEGs of cold atomic gases offer good opportunities to 
explore novel properties of Bogoliubov excitations. 



Recently, Kovrizhin and co-workers 



m 



10| predicted that the Bogoliubov mode exhibits 



striking tunneling properties. They showed that the transmission probability of Bogoliubov 
phonon through a potential barrier increases in the low energy region with decreasing the 
incident energy. In the low-energy limit, the perfect transmission is realized irrespective of 

Bogoliubov mode is referred 



the height of the barrier. This interesting tunneling property o 
to as the anomalous tunneling Since their prediction 9, 



10| , the anomalous tunnelin, 



has attracted much attention, and has been addressed by many papers 
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16 



171, [181, Hi . 



13|, 



14, hi 



As the origin of the anomalous tunneling effect, various mechanisms have been proposed, 
such as quasiresonance scattering 1^ , localized components of Bogoliubov mode appearing 



near the barrier 



ll|, and anomalous enhancement of quasiparticle current [12]. For the 



perfect transmission in the low-energy 
densate and excitation wave functions 



Bogoliubov phonons llJl, has been pointed out. The anomalous tunneling phenomenon was 
shown to occur even in the supercurrent state [HI, uM, 12Q| , as well as at finite temperatures 
It has been also studied in the presence of a periodic potential [l5|, ilGj, as well as a 



imit, the importance of the coincidence of the con- 
as well as supercurrent behavior of low-energy 



random potential [17[. It has been also pointed out that similar phenomena to this can be 
seen in the scattering of Bogoliubov phonon by a spherical potential in three dimensions 
djisl- as well as the refraction of Bogoliubov phonons lol. 



In this paper, we investigate tunneling properties of Bogoliubov phonon in a BEG at 
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T = 0. In Ref. [iz], we have considered the case when the incident and transmitted Bogoh- 
ubov phonons feel the same condensate densities on both the right and left of the barrier. 
In this paper, we extend this previous paper to the case when the condensate density is 
different between the right and the left of the barrier. As an extreme case, we also deal with 
the case when the condensate density is absent on one side of the barrier. Applying the 
finite element method to the Bogoliubov coupled equations, we numerically calculate the 
transmission probability and phase shift of Bogoliubov phonons. We find that Bogoliubov 
phonons twist the phase of the BEC order parameter (condensate wave function) due to a 
back-reaction effect, which leads to the induction of Josephson supercurrent. The induced 
supercurrent is shown to satisfy the Josephson relation with respect to the twisted phase 
when the condensate density is the same on both sides of the barrier. The supercurrent 
is induced only in the region near the barrier when the condensate has the same densities 
across the potential barrier. In the case when the condensate density is different between 
the right and left of the barrier, the supercurrent is also induced in the region far away from 
the barrier. In addition, the excess quasiparticle current is supplied from the condensate to 
conserve the total current, so that one obtains the enhancement of the transmission proba- 
bility of quasiparticle current in the low-energy region. We also show that the supercurrent 
is induced when one injects free atoms from the outside of condensate. In addition, atoms 
are shown to evaporate from the surface of superfluid region when Bogoliubov excitation 
propagates toward the superfluid-normal phase boundary. 

This paper is organized as follows: in Sec. [Tll we present the model and formalism of the 
Bogoliubov mean-field approximation, as well as the finite element method which we apply 
for solving the Bogoliubov equations. In Sec. IIIIl we study the tunneling of Bogoliubov 
phonons through a rectangular potential barrier. We give a detailed discussion on the origin 
of the anomalous tunneling and induced Josephson supercurrent. In Sec. HVl we study the 
tunneling in the presence of a step potential which yields the different condensate densities 
across the potential barrier. In Sec. |Vl we discuss the tunneling of excitations and atoms 
between superfluid and normal regions. 
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FIG. 1: Schematic of the system. The arrows on the left describe the incoming (upper arrow) and 
reflected (lower arrow) Bogoliubov excitations. The arrow on the right describes the transmitted 
one. 

II. MODEL AND FORMALISM 



We consider tunneling phenomena of Bogoliubov mode through a barrier potential, as 
schematically shown in Fig. 1. We assume that the barrier potential U{x) only depends 
on X and ignore the motion of atoms in the y and z directions, so that we consider a one- 
dimensional tunneling problem along the x direction. This kind of one-dimensional geometry 
has been recently realized 21\. In Ref. |2ll], a BEC was prepared in a narrow elongated trap 
with a wall-type potential barrier, which varies only in the axial direction and the potential 
width is much longer than the radial size of the gas cloud. We also ignore temperature 
effects as well as effects of a harmonic trap. The latter assumption is justified when the 



BEC is trapped in an elongated trap 



22| or a box-shaped trap [23 1 



We treat the tunneling of Bogoliubov mode within the Bogoliubov mean-field theory for 



a weakly interacting Bose gas at T = 12: 



251]. To describe the BEC phase, we divide the 



Bose field operator '?/'(r) into the condensate wave function ^E'ol^) and the noncondensate 
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part, as 



j 



(1) 



where a\ is the creation operator of a Bogohubov excitation in the jth state. The condensate 



wave function "^oi^) = satisfies the static Gross-Pitaevskii (GP) equation 



24 



(2) 



Here, m, /i, and U{r) represent the mass of a boson, chemical potential, and barrier potential, 
respectively. g{> 0) is a repulsive interaction between bosons. In Eq. ([1]), Uj{r) and Vj{r) 
satisfy the Bogoliubov coupled equations. 



-^ + U{r) + 2g\^o\'-f^ 
-^ + U{r)+2g\^o\''f^ 



Uj - g'^lvj = EjU 



'31 



V, - g{%)\ = -E,v 



J' 



(3) 
(4) 



where Ej is the Bogoliubov excitation spectrum. To solve Eqs 



and (jlj) with an appro- 



priate boundary condition, we use the finite element method [27|. For this purpose, it is 
convenient to rewrite Eqs. ([3]) and (jlj) in the matrix form 



where 



H 



H - E,r, 



-V2 + f/(r) + 2|^oP-l 



0, 



*\2 



-V2 + f/(f) + 2|*op-l 



^3 




(5) 



(6) 



(7) 



(8) 



In Eq. ([5]), we have introduced dimensionless variables r = r/^, Ej = Ej/fi, U = U/fi, 
and ^'o = ^o/'^o? where no = fi/g is the condensate density far away from the barrier, and 



= 1/ y/2mgnQ is the healing length. To simplify our notations, we omit the bars and indices 
of eigenstates in the following part of this section. Equation ([5]) can be obtained from the 
variational principle 6L = 0, when the Lagrangian L has the form 



L = [ dr [(V0t) • (V0) + 0^ (f/(r) + 2|^or - 1 - ^o^+ - (^S)'^-) - Er^] 



(9) 



Here, Q is the volume of the system, and t± are given by 





(10) 



We introduce spatial positions r-j (i = 1, 2, . . . , N\jn the system, which are referred to 
as nodes in the literature of the finite element method {2?]. We then assign the interpolation 
function Ni{r) at each r,, which equals unity at r = rj and linearly decreases to zero at 
adjacent nodes of rj. Namely, the interpolation function satisfies 



For example, in one-dimensional case, we define the nodes at Xi {i 
interpolation function Ni{x) is given by 



1,2, 



^z+1 



tjC 2 ]^ iC tjC 1} 



(Xi < X < x^ 



(11) 

AT). The 



(12) 



0, {x ^ Xi^X ^ Xj+i). 

Using Ni{r), one can approximately write u{r) and v{r) in the forms 

u{r) = ^UiNi{r), 

i 
i 

Substituting Eqs. (IT^ and ffUl) into Eq. Q, we obtain 



(13) 
(14) 



(15) 
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where 



(16) 



K,,, = I dr (ViV,) ■ (ViV,) , 

I dr Ni [E ± (t/(r) + 2|^o|' - 1)] Nj 
(EtI) I dr NiN, ± f/^ /" dr N^N^Ni 
±2 J2ifi - m){fi' + ^91') I dr N.N.NiNv, (17) 



Pij = / dr Ni^lNj 



= Y^^fi + Wi){fi' + Wi') I dr N,N,NiNi,. (18) 

In obtaining Eqs. f|T7|l and (|T8|l . we have expanded U{r) and '^o{r) as 

U{r) = J^^iNiir), (19) 
Mr) = Y,Ui + Wi)Ni{r). (20) 

Here, f/, = f/(n), = Re[^o(n)], and gi = Im[^o(n)]. In Eqs. dH-dUD, fi^j, fii,,,;, and 
mean that the integrations are carried out in the regions where NiNj, NiNjNi, and 
NiNjNiNii are finite, respectively. The integratioiis in Eqs. f|T6|) - f|T8|) can be evaluated in 



the standard manner of the finite element method 



Equation ([T^ can be rewritten in the matrix form as 

L = u^ (^k- M") u + v^ (^k + M+) V - uPv - v^P*u, (21) 

where {u}i = Ui, {v}i = Vi, {k}i j = Kij, {M'^}ij = M^j, and {P}i,j = Pij- The equations 
for u and v are, respectively, obtained from 5L/5u^ = and 5L/5v'^ = 0, which give 



K - EM'j u-Pv = 0, (22) 
k + EmA v-Pu = 0. (23) 



The advantage of using the finite element method is that one can obtain the solutions by 
simply diagonalizing Eqs. (12^ and (12^ under an appropriate boundary condition, instead 
of solving the differential Eqs. ([3]) and (jlj). In the following sections, we will numerically 
solve Eqs. (I22l) and (l23l) for given barrier potentials. 



III. TUNNELING THROUGH THE RECTANGULAR POTENTIAL BARRIER 



In this section, we consider the one-dimensional tunnehng problem of Bogoliubov excita- 
tions through a rectangular barrier potential shown in Fig. [H The potential barrier is given 
by 

U{x) = Uo9(^^- \x\^ , (24) 

where 6{x) is the step function. Uq and d describe the height and width of the barrier, 
respectively, and we consider the case of repulsive potential barrier {Uq > 0). In this section, 
we treat the case when the condensate densities are the same on both sides of the barrier, 
as shown in Fig. [H Although this case has been examined in our previous paper 12] , we 
give further analyses for the tunneling of Bogoliubov phonon here. In Sees. |IV] and |Vl we 
will also compare the results in this section with the case when the condensate density on 
the left of the barrier is different from that on the right of the barrier. 
In the present case, the GP equation can be solved analytically [l^, as 



tanh 



1 



x\ 



+ arctanh7 



x,q 



(|x| >d/2), 
{\x\<d/2). 



(25) 



Here, d = d/^, K = (5'^ + 2(f/n - 1), and q = K/y^K^TW, where Uq = Uo/fi. (cn(x,g) 
is the Jacobi's elliptic function [28].) 7 = ^q{x = d/2) and f3 = 'I'o(O) are determined from 
the boundary conditions in terms of \l/o(a;) and d'^^/dx at x = ±(i/2, which give 

/5 



7 



cn 
1 



2 2' " 



7 



2f/n 



(/3^ + 2(t/o-l)/52 + l). 



(26) 
(27) 



The values (3 and 7 are determined by numerically solving Eqs. (!26|) and (127|) . 

To solve Bogoliubov equations fl22l) and fl23l) . we need asymptotic solutions for x = ±00. 
In our tunneling problem, each eigenstate with index j in Eqs. ([3]) and (jlj) corresponds to 
Bogoliubov excitation with energy E injected from one side of the barrier. In Sees. IIIHTV] we 
omit the index for eigenstates for simplicity. Far from the barrier (|x| C,), the Bogoliubov 
mode is described by the plane-wave {u{x),v{x)) = {uE,VE)e^P^. Substituting this into 
Eqs. (|3]) and (jlj), one obtains the well-known Bogoliubov excitation spectrum as l| 



E„ 



(28) 



where Sp = /2m. Namely, for a given mode energy there are four particular solutions 
in terms of the momentum p, given by 



P 



±\p2im\l ^/WT~(gnof - gno = ±k, 



±iy/2m\/ \/ ^ {gnoy + gno = ±iK. 



(29) 



The first two solutions {p = ±k) describe the ordinary propagating waves in the in- 
directions. The remaining two imaginary solutions {p = ±m) describe localized states. 
We note that while the latter localized solutions are actually not necessary in a homoge- 
neous system, we cannot ignore them in the present inhomogeneous system. The amplitudes 
of the propagating components are given by 

a 




+ 1 



E 



(30) 



where L is the system size in the x direction. On the other hand, the amplitudes for 
the localized states are given by {u\,v^) = (—6, a). Thus, in contrast to the propagating 
solution in Eq. fl5Ul) . the normalization of the localized components becomes negative as 

(u'^Er - K)' = -yL. 

Using the propagating solution (m^, t>^)e^*'^^ and localized one {u^, T;^)e^''^, we construct 
the asymptotic forms of the Bogoliubov wave function for x ±oo. Assuming that the 
Bogoliubov phonon is injected from x = — oo, we obtain the asymptotic solutions as 





(u\ 






< 


)"< 




u 








( \ 

a 



vV 



+ A 



Akx 



B 



1-^ 



ya J 



[x — oo 



(31) 



oo) 



Here, r and t are, respectively, the reflection and transmission amplitudes, which satisfy 

1. (32) 



As will be discussed later, this condition is deeply related to the conservation of energy 
flux. In Eq. (13T|) . A and B represent the amplitudes of the localized components near the 
potential barrier. 
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FIG. 2: Calculated transmission probability W and phase shift 6 as functions of the incident 
energy E for a rectangular potential barrier. We set the width d and height Uq of the barrier as 
{d,Uo) = (C, 2/i) (solid line), (CjS//) (dashed line), {^,10fi) (dotted line), and (4^,2/i) (dash-dotted 
line) . 

We numerically solve the Bogoliubov coupled Eqs. ( l22l) and ( l23l) for a given incident 
energy E. In this procedure, the condensate wave function in Eq. (123]) is used, and the 
solution is determined so as to satisfy the asymptotic solution in Eq. (pTl) . 

Figure [2] shows the calculated transmission probability W = as well as phase shift 
S = arg(t), as functions of the incident energy E for various barrier heights and widths. We 
call attention to the characteristic features of W and 6 in the low-energy region (E/ n < 0.5). 
One can clearly see the anomalous tunneling behavior discussed in js, 9, 10 1 in Fig. [21 



Namely, below a certain incident energy {E/fi ~ 0.5), W increases and 6 decreases with 
decreasing E, in contrast to the behaviors above that energy [W decreases and 6 increases 
as E decreases). Furthermore, W and S approach unity and zero in the low-energy limit 
E —>■ 0, respectively, irrespective of the values of d and Uq. When the incident energy E is 
very large {E ^ fi), since the Bogoliubov phonon loses its collective nature, the tunneling 
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FIG. 3: Spatial variation in the Bogoliubov wave function (u(x),f(x)). (a) E/^ = 0.01 ^ 1 
(anomalous tunneling). We set (d,C/o) = (^,10^). (b) E/ix = 1.68. We set ((i,C/o) = (4^,2^). In 
this case, the resonance tunneling {W = 1) is realized, as shown in Fig. [2] (see the dash-dotted 
line). The dotted line is the potential barrier U{x) in units of /x. 

property becomes close to that of a single particle. 

We note that the perfect transmission of Bogoliubov phonon {W ^ 1,5 — > 0) shown in 
the low-energy limit in Fig. [2] is quite different from the typical tunneling properties of a 
single particle, where W and 5 approach and — 7r/2 in the low-energy limit, respectively. 
Namely, in the latter case the particle is completely reflected by the potential barrier {29! . 

We also note that the energy region in which W and 5 exhibit the anomalous tunneling 
behavior {W increases and 5 decreases with decreasing E) depends on the height Uq and 
width d of the potential barrier. This region becomes narrower for higher and wider potential 
barrier, as shown in Fig. [2l 

In Fig. [2], we find that 11^ = 1 is also obtained at finite energy {E/ ^ = 1.68) in the 
case of {d,Uo) = (4^,2yu), due to the resonance tunneling effect. To see the difference 
between the resonance tunneling effect and the anomalous tunneling effect, we show in 
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Fig. [3] the wave functions in the two cases. In the case of resonance tunnehng, one sees that 
while |up is enhanced in the barrier, |f p is suppressed there. The peak of \u\^ is a clear 
signature of the formation of a resonance state. The suppression of |f p indicates that the 
Bogoliubov excitation behaves like a single particle during the tunneling through the barrier. 
In contrast, in the case of the anomalous tunneling, both \u\'^ and |f p simply become small 
in the barrier and almost coincide with each other. Indeed, it was shown in Refs. Q, [l^ that 
u{x) and v{x) reduce to the condensate wave function '^q{x) in the low-energy limit. The 
difference mentioned above indicates that the anomalous tunneling and resonance tunneling 
are different phenomena. 

We briefly note that, as shown in Ref. the anomalous tunneling effect originates 
from the fact that the wave functions of a Bogoliubov phonon with a small momentum p 
has the same form as the condensate wave function in the supercurrent state, accompanied 
by a finite superfiow Jg = nop/m. Recently, Morgan et al. [30] have presented a modified 
Bogoliubov theory where the wave function of Bogoliubov mode is constructed so as to be 
orthogonal to the solution obtained from the GP equation. Since their formalism does not 
affect the current-carrying component of the Bogoliubov wave function (which dominates 
the anomalous tunneling phenomenon), the perfect transmission of low-energy Bogoliubov 
phonon is still expected to occur. Thus, the anomalous tunneling phenomenon does not 
depend on the definition of the wave function of Bogoliubov mode. 

Propagation of Bogoliubov phonon is accompanied by quasiparticle current Jq, as well as 
energy flux Qq. When one uses the asymptotic solutions in Eq. ( l3Ti) . they are given by 



^(1 - |r|2) - l^e'^^ {k Im [A{e-'^'' + r*e'^'')\ + Re [A{e-'^'' - r*e'^'=)\) , 

{x « -0, (33) 

- ^e"''^ {k Im [tB*e'^''] +kRe [tB*e'^'']) , (x > 0, 

g^H"'"^"'^'"""''- (34) 
!§{a' + l?)\t\\ 

The detailed definitions of Jq and Qq are summarized in Appendix B. Since the energy 
flux Qq is conserved (see Appendix B), one obtains |rp + |tp = 1 from Eq. f lMl) . Using 
this, we find that the quasiparticle current Jq is also conserved in both limits x = ±oo as 
Jq(x = — oo) = k{l — \r\'^)/mL = Jq(x = oo) = k\t\'^/mL. However, except for the limits 
X = ±oo, the last terms in Eq. (|33|1 become finite, which come from the coupling between 

12 



the propagating and localized components in Eq. (l3Ti) . As a result, while Qq is conserved 
everywhere, we expect that Jq is not conserved near the barrier. 

To see the non-conserving behavior of Jq, we directly evaluate it using the solution of 
Bogoliubov equations fl22]) and fl25]) . As shown in Fig.|H^a), we obtain the excess quasiparticle 
current 

AJq(x) = Jq(x) — Jq(x = — Oo) (35) 

near the barrier. Namely, when the Bogoliubov phonon approaches the barrier, Jq is en- 
hanced. Jq is constant in the barrier, and it decreases to be Jq(oo) = Jq(— oo) when the 
phonon goes away from the barrier. In Fig. |4](a), the enhancement of AJq occurs near the 
barrier where the condensate density ns{x) deviates from no[= ns{x = ±00)]. 

In Fig. 111(a), we find that the enhancement of AJq is more pronounced for lower incident 
energy E. In addition, as shown in Fig.[5l the excess quasiparticle current is more remarkable 
when the barrier is high, although the energy region where AJq(x = 0) is large is narrower 
for larger Uq. Since AJq(x = 0)/{k/mL) approaches a constant value, we also find from 
Fig. [5]that AJq(x = 0) is proportional to the incident momentum k in the low-energy limit. 
[Note that AJq(x = 0) in Fig. [5] is normalized by the incident quasiparticle current k/mL] 

The enhancement of quasiparticle current near the potential barrier implies that more 
quasiparticles than those carried in the incident current impinge on the barrier. Apparently, 
this is expected to lead to the increase in the transmission probability of quasiparticles. 
Indeed, comparing the result for (d, Uq) = (^, 10/i) in Fig. H] with the corresponding result in 
Fig. [21 one finds that the energy region where the anomalous enhancement of transmission 
probability is obtained {E < O.lyu) coincides with the region where the excess quasiparticle 
current AJq(x = 0) is remarkable. 



As shown in Ref. [12|, the excess quasiparticle current is supplied from the condensate. 
Namely, the transmission of Bogoliubov phonon is considered to be assisted by the supply 
of excess current from the condensate. Thus, in a sense, the mechanism of the anomalous 
tunneling may be considered as a kind of screening effect by Bose condensate. This argument 
partially explains the physical mechanism of the anomalous tunneling effect discussed in 
Refs. jsiSllol. However, apart from the enhancement of low-energy transmission probability, 
this argument is not enough to explain the perfect transmission in the low-energy limit. In 
this regard, in Ref. 1^, we have shown that the perfect transmission can be understood as 



a result of the supercurrent behavior of low-energy Bogoliubov phonon. 
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FIG. 4: (a): Excess quasiparticle current AJq{x) = Jq{x) — Jq(x = — oo) when (d^Uo) = 10^). 
The dash-double dotted hne is the condensate density ns{x) in units of uq. (b): Phase 0{x) 
of the condensate wave function ^'0(2;) created by the tunnehng of Bogohubov phonon. We set 
9{x = —00) =0. (p is the phase difference between condensates at x = ±00. Nq = uqL is the 
number of condensate atoms. In both the panels, the dash-dotted hne indicates the region of the 
potential barrier. 



In Ref . 12] , it was found that the counterflow of supercurrent is induced near the potential 
barrier due to a back-reaction effect of quasiparticle current, which restores the conservation 
of total current. The induction of supercurrent indicates that the phase of the BEC order 
parameter \l/o(a;) is twisted by quasiparticle current as 



(36) 



(Here, we assume that the amplitude of the condensate wave function is unchanged.) The 
induced supercurrent by this phase modulation is given by 



AJJx) 



-dj(x) 



m 



(37) 
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FIG. 5: Excess quasiparticle current in the potential barrier AJq(x = 0) as a function of the 
incident energy E. 



As shown in Appendix B, AJs{x) is related to the excess quasiparticle current AJq as 



AJJx) 



-AJq(a;) 



(38) 



[Here, we set (ajoj) = 1 in Eq. (lB3ip assuming that one Bogoliubov excitation is injected.] 
As a result, the phase 9{x) is evaluated to be 

6{x) = —m 



dx — - 



(39) 



Namely, the phase modulation is caused by the excess quasiparticle current AJq. The as- 
sumption in Eq. (l36ll is valid as long as 6{x) is small, because the change in the amplitude 
of the condensate wave function gives higher-order corrections. Since 9{x) is inversely pro- 
portional to the number of condensate atoms Nq, as shown below, 9{x) is negligibly small, 
so that the assumption in Eq. (!36|) is justified. As discussed in Appendix [Bl the inclusion of 
the back-reaction effect of quasiparticles on condensates requires the modification of the GP 
equation as Eq. (1B29I) . In the present case, the new condensate wave function including the 
back-reaction effect is perturbatively obtained with use of the ansatz in Eq. ( l36i) without 
solving Eq. flB29l) . 
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FIG. 6: Induced supercurrent AJs{x = 0) as a function of the relative phase (p across the potential 
barrier. The slopes of the lines are 7 = 0.14, 0.050, and 0.015 for {d,Uo) = (^,5/i), and 

(^,10/i), respectively. 

Figure ID^b) shows 6{x) when the Bogoliubov phonon is injected from x = —00. The 
spatial variation in the phase 6{x) is remarkable near and in the barrier, where large excess 
current AJq is obtained [see Fig. |l](a)]. 

Figure [6] shows the magnitude of induced supercurrent at x = as a function of = 



^(-00) - e{x) [see Fig. SKb) 



Josephson current relation [311 



We clearly see that |AJs(x 



= I J sin ~ Ijc, 



«1). 



satisfies the ordinary 



(40) 



(In our case, since is proportional to the inverse of the total number of Bose-condensed 
particles A^^o, so that 0^1.) The Josephson critical current Ij in the present case has the 
form, 



no 
mi 



(41) 



where 7 is determined from the slope of the lines in Fig. [6l This result means that the 
Josephson critical current Ij may be evaluated from the analysis of quasiparticle tunneling 
without directly examining the Josephson current. 
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Finally, we remark that the tunneling properties of Bogoliubov excitation discussed in 
this section suggest an important role of Bogoliubov phonons on the fluctuation of the 
relative phase between two condensates at finite temperatures. When Bogoliubov phonons 
are excited on both sides of the barrier at finite temperatures, they tunnel through the 
potential barrier and twist the relative phase. This is expected to lead the fluctuation of the 
phase difference between the condensates on the left and right of the barrier. In particular, 
large phase fluctuations may be induced in the temperature region where the population 
of Bogoliubov phonon becomes dominant. This phase fluctuation due to the tunneling of 
Bogoliubov phonons could be observed in a BEC in a double-well potential, where the 
thermally induced fluctuations of the relative phase between two condensates were recently 



observed 



32|. 



IV. TUNNELING BETWEEN CONDENSATES WITH DIFFERENT CONDEN- 
SATE DENSITIES 

In Sec. Illlt we considered tunneling properties of Bogoliubov phonons through the rect- 
angular potential barrier in the case when the left and right of the barrier have the same 
condensate densities. In this section, we consider the more general case when the condensate 
densities are different between the right and left of the barrier. This situation is achieved 
by simply imposing a uniform potential on the right side of the barrier, as 

U{x) = Uoe (^^ - \x\^ + U,e (a: - . (42) 
In this case, the condensate density at x ^ oo is given by 

~ T . ^2 (0<f/i</i), 

no = ^o(x = ^f= I' (43) 

[O, (f/i>/i). 

In this section, we consider the case of < f/i < /i. The case of f/i > /i will be discussed in 
Sec.lVl The barrier potential, as well as the condensate wave function '^q{x), is schematically 
shown in Fig. [71 

To solve Bogoliubov equations ( l22i) and (!23ll . we construct the condensate wave function 
\l/o(x), as well as the asymptotic solutions at x = ±oo. The former is analytically obtained 
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u, 







U(x) 



u 



X 



FIG. 7: Schematic of the system in the presence of the barrier potential in Eq. ([42 



from Eq. ([2]) as 



tanh 



V2 



4^ ( a; + ^ ) + arctanh72. 



{x-X(i),q 



\/l — Ui tanh 



(x - f ) + arctanh ( 



{\x\ < d/2) 
, (a:>f), 



where ?7i = 7l = ^o(— c^/S), and 7ij = \l/o(c?/2). xq satisfies the conditions \I/(xo) 

and (i\l/o(a;)/(ix|a;=i,.(, = 0. Xq, /3, 7^,, and 7^ are determined from the equations, 

IR = - 



cm ,/^(|-xo),g 



7l 



7I 



7i 



/9 



cm ,/^(| + xo),g 



1 



(/?^ + 2([/o-l)/32 + (l-t/i)^), 



2(f/o - f/i 

(/?^ + 2(t/o - l)(3' + 1) . 



Equations (H5]) - fl48p are derived from the boundary conditions at x = ±d/2. 
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The asymptotic solutions of the Bogohubov equations at x = ±00 are obtained in the 
same manner as in Sec. IIIII Assuming that the Bogohubov phonon with the energy E = 
\/ep{ep + 2gh) is injected from x = —00, we have 





/ \ 
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J 




u 




[vj 



1 \ 

a 






w 


v 



/ \ 



a 
b 



-ikx 



+ A 



ya J 



-KnX 



[X 



-00 



(49) 



{x — i> 00). 



Here, k, k, and (a, 6) are given in Eqs. (l29l) and ( l30l) . The parameters appearing in the 
asymptotic solution at x = 00 are given by 



kR 

I^R 




2L 



+ 1 



v 



2L I E 



\ 



J 



(50) 
(51) 

(52) 



Using the condensate wave function in Eq. (^^, we numerically solve the Bogohubov 
Eqs. ( 122|) and ( l23i) so as to satisfy Eq. ( H9i) . Once the wave function {u{x),v{x)), as well as 
r and t, are determined, we can calculate the transmission probability from the conserving 
energy flux Qq. The energy flux Qq at x ^ ^ is given by 

kRE 



mL 



(53) 



The transmission (reflection) probability W (R) is conveniently defined as the ratio of the 
incident and transmitted (reflected) components of Qq. From Eqs. (IMI) and (153|) . we obtain 



W 
R 



kRjaR + hi) 2 
k{a? + 62) II' 

l„|2 



(54) 
(55) 



Equations flMl) and flS5]) satisfy the relation R + W = 1 because of the conservation of Qq 
as proved in Appendix [B] 

We note that, when we calculate the transmission probability from the quasiparticle 
current, we obtain a different result from Eqs. (|5^ and (!55|l . Using the expression for the 
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quasiparticle current at x ^ ^, 

k In h 

Jq = ^e-"«"(/€^jIm [tB*e'^^''] + fc^Re [tE*e'^«"]), (56) 

iTtij in 

and Eq. (l33l) . we define the "transmission (reflection) probability" VTj (-Rj) as the ratio of 
the incident and transmitted (reflected) components of Jq(x = ±00). Then, we flnd 

Rj = |r|2 = R. (58) 

Since Wj > W, Eqs. ([57j) and (jSHD do not satisfy the condition Rj + Wj = 1, unless f/i = 0. 
This is because of the fact that Jq is not conserved, as discussed in Sec. [Till and Appendix 
[B1 When Ui = (this case was discussed in Sec. IIIII) . the breakdown of the conservation of 
Jq is restricted to the region near the barrier. Namely, all the supplied component AJq is 
completely absorbed after the quasiparticle is transmitted in the right condensate, as shown 
in Fig.m As a result, the transmission probability, which is deflned using Jq(x = ±00), is not 
affected by this non-conserving character of Jq. On the other hand, the fact of Wj + Rj > 
when Ui > indicates that the non-conserving behavior of Jq remains even at a; = cxo. 

Figure [8] shows the calculated transmission probability W, Wj in Eq. fl571) . and the phase 
shift 6 = arg(t). While the phase shift 6 approaches in the low-energy limit irrespective 
of the value of f/i, the perfect transmission {W — * 1 in the low-energy limit) is absent when 
Ui > 0. In Ref. [191], Watabe and Kato have obtained the analytic expressions for W and S 
in the low-energy limit for arbitrary potential barrier shape. According to their results, W 
and 6 read [lol 

4v^r^ ^^ 5^0, (E^O). (59) 

(1 + ^/1^)2 



These results can be also obtained in the case of a (5-function potential barrier [33|]. Our 
results in Fig. [8] are consistent with their earlier results in Eq. fl59|) . Equation fl59l) shows 
that the transmission probability W becomes less than unity when Ui > 0. As pointed 



out in Ref. [19[ , the low-energy behaviors of W and 6 are determined only by the potential 
difference at a; = ±00 {Ui in the present case), and they do not depend on the detail of the 
potential barrier in the middle. 

In Fig. [8](b), one flnds that Wj is remarkably enhanced to be larger than unity in the 
low-energy region. To see the relation of this large enhancement and the non- conserving 
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FIG. 8: Transmission probability W (a) calculated from the energy flux, and Wj (b) calculated 
from the quasiparticle current Jq, and phase shift 6 (c). We set (d, C/q) = 5/x). 

character of Jq, we show the spatial variation in Jq in Fig. [9](a). Comparing this result 
with Fig. m we find that the excess quasiparticle current A Jq remains finite even far away 
from the barrier (x ^ d) when Ui > 0. This excess current is found to be supplied from the 
condensate through the source term S{x) defined in Eq. (lB14p . as shown in Fig. [n](b). Figure 
Mjo) also shows that this supply dominantly occurs in front of the barrier (—5 < x/^ < 0). 
(Note that the phonon is injected from x = —oo.) As a result, Wj given by the ratio of 
incident and transmitted quasiparticle current is remarkably enhanced. 

As discussed in Sec. Illlt the excess component A Jq(x) = Jq(a;) — Jq(x = — oo) is cancelled 
out by the counter flow of supercurrent to conserve the total current. As shown in Fig. M^c), 
the induced supercurrent remains finite even at a; ^ oo, which is in contrast to the case 
of Ui = 0, where AJg is only finite near the barrier. The reason for this can be considered 
as follows: as discussed in Sec. Illlt Bogoliubov phonons twist the condensate phase when 
they tunnel through a potential barrier. In addition, Bogoliubov phonons can be regarded 
as quantized oscillations of the phase of the condensate wave function j^. Since the phase 
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FIG. 9: Spatial variation in quasiparticle current Jq{x) (a), source term S{x) (b), and induced 
supercurrent AJg (c) when E = 0.01/i. We set {d, Uq) = 5/i). The dash-dotted hne in (a) shows 
the potential barrier U{x) in units of /i when Ui = 0.5fi. 

stiffness is weak on the right side of the barrier due to the small condensate density, the 
transmitted Bogoliubov phonons can easily twist the phase of the right condensate when 
Ui > 0. This leads to the induction of counter superfiow far away from the barrier. Indeed, 
\AJs{x ^ ^)| and Jq{x ^ ^) are larger for larger Ui, as shown in Fig. [9l 



V. TUNNELING BETWEEN SUPERFLUID AND NORMAL REGIONS 

In this section, we consider the case when the condensate density at x ^ — ^ is absent. 
To realize this situation in a simple manner, we use the potential 

U{x) = U2e{-x) (60) 

with U2 > /i. In what follows, we call the negative x side the normal region and the positive 
X side the superfiuid region. In the normal region, Bogoliubov excitations reduce to free 



22 



atoms, having the energy 

E; = ep + {U2~ (61) 

Here, we discuss two different tunnehng problems, i.e., tunnehng of atoms from the normal 
region to the superfluid region (N-S tunneling), and the tunneling of Bogoliubov excitations 
from the superfluid region to the normal region (S-N tunneling). These two cases enable 
us to study how free atoms are injected into a condensate and emitted from the surface 
of the condensate. We note that these tunneling problems are analogous to the quantum 
evaporation and condensation at a free surface in superfluid ^He 34 1. 

We first consider the N-S tunneling. The analytic solution of the GP Eq. (|2]) is given by 

A 



sinh( A^-C 



(x < 0), 



tanh ( + arctanha ) , (x > 0) , 



(62) 



where a = l/v^, A = ^2{U2 - 1), and C = log[(A + Va^ + A^)/^], 
The asymptotic solution of the Bogoliubov equations has the form 





( \ 
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< 






u 




[vj 




[X 



-oo 



(63) 



oo 



where k, k, and (a, 6) are given in Eqs. fl2^ and fl5U]) . The wave numbers and kl for 
X ^ —oo are given by 

kL = V 

kl = V2m^/E + (f/2 - fx). 



(64) 
(65) 



kl and kl are propagating and localized waves for x — oo, which are obtained by solving 
E = ±Ep in terms of p. We note that the localized v component in Eq. ( I63i) describes the 
(proximity) effect of the condensate in the normal region. 

Using Eq. (1631) . we obtain the quasiparticle current Jq, as well as the energy flux of 
quasiparticles Qq, in the normal region {x —^) as 



^ mL 
kiE 

mL 
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(i-kH- 



(66) 
(67) 





k , 


Wj = 




Rj = 


|r|l 



Since the localized v component in Eq. (163|) does not give rise to any contribution to the 
currents, Jq and Qq reduce to those of free atoms which satisfy the relation Qq = EJq. From 
Eqs. and (p7|) . we obtain the transmission (reflection) probability W (R) as 

W = L-^{a^ + b^)\t\^, (68) 
R = |rp. (69) 

Since the energy flux is conserved as shown in Appendix [B], they satisfy the condition R + 
W = 1. We also obtain the "transmission (reflection) probability" Wj {Rj) for quasiparticle 
current from Eqs. (!33l) and (!66l) . 

(71) 

We again find that, W and Wj do not coincide with each other. Because of Wj < W, we 
obtain Rj + Wj < 1. This implies that the quasiparticle current decreases on the superfluid 
region. Since the ratio between W and Wj is given 
and Wj become equal when E/gno ^ 1. 

Figure [TO] shows the transmission probability W obtained from the energy flux and the 
phase shift 6 = arg(i)- We also show the quasiparticle transmission probability Wf^ in 
Fig. [11] (a). In Figs. [TOjandfTTl we note that the origin of E is taken to be U2 — /i, because 
atoms are perfectly reflected when E < U2 — fJ^, leading to vanishing Wf^ and W. We 
find that both W and Wj decrease with decreasing E, and W and Wj approach when 
E ^ U2 — fi, while 6 approaches a positive value in the limit of E ^ U2 — fi- Thus, the 
anomalous tunneling behavior does not occur in the present case. 

Figure fTWa) shows the spatial variation in quasiparticle current Jq{x), source term S{x) 
[defined by Eq. flB14l) ]. as well as the induced supercurrent AJs(a;). The existence of trans- 
mitted component of Jq(x) shows that the incident current of free atoms from the normal 
region is converted into the Bogoliubov excitations inside the condensate. Furthermore, one 
finds that Jq{x) decreases near the surface at x ~ 0, and the supercurrent AJs{x) is induced 
around the surface. The source term S{x) becomes negative near the surface of the superfluid 
region reflecting the behaviors of Jq{x) and AJs{x). These phenomena indicate that injected 
atoms are Bose condensed in the superfluid region, which give rise to the supercurrent A Jg. 
The condition Rj + Wj < 1 reflects the fact that a part of the incident current of free atoms 
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FIG. 10: Transmission probability W obtained from the energy flux (a) and phase shift S (b) 
both the N-S and S-N tunnehng cases. 
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FIG. 11: Transmission probabihties obtained from the quasiparticle current in the N-S {Wf ) ( 

and S-N (W^^) (b) tunnehng cases. 
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FIG. 12: Spatial variation in quasiparticle current Jq{x) (solid line), induced supercurrent AJs{x) 
(dashed line), and source term S{x) (dotted line), when U2 = 2/i and E = l.Ol/U for the N-S 
tunneling (a) and S-N tunneling (b). Current density and source term are in units of ki/mL and 
kL/mL^ in (a), and k/mL and k/mL^ in (b), respectively. The dash-dotted line in (a) indicates 
the condensate wave function '^q{x) in units of -^/rio/0.3. 



is converted to supercurrent inside the condensate. The supercurrent AJs decreases as E 
increases, because the character of produced Bogohubov phonon becomes close to that of 
single-particle excitation, as E increases. As a result, Wj approaches W when E ^ quq. 

We next consider the S-N tunneling. Assuming that the incident Bogoliubov mode comes 
from X = +00, the asymptotic solutions {u,v) for x —>■ ±00 are given by 
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(72) 



00 



Using Eq. ( 1721) . the quasiparticle current Jq and energy flux of quasiparticles Qq in the 
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normal region (x ^ — ^) are calculated as 

J, ^ ^If , (73) 

From Eqs. and ([71]), we obtain the transmission (reflection) probability W {R) obtained 
from the energy flux as 

R = |r|2. (76) 

Equations (1751) and (1751) again satisfy the condition i? + = 1. From Eqs. (15^ and (1751) . 
we obtain the transmission (reflection) probability Wj {Rj) for quasiparticle current as 

Wj = '^\t\^ = L{a^ + b^)W, (77) 
Rj = \r\'. (78) 

In contrast to the N-S tunneling, it is clear from Eq. (1771) that Wj > W, which leads to the 
condition Rj + Wj > 1. This implies that the quasiparticle current is supplied around the 
surface at x ~ 0. When E ^ guQ, Wj reduces to W. 

It can be generally shown for the Bogoliubov coupled Eqs. ([3]) and (jlj) that W and 5 are 

ler the incident wave comes from x = —oo (N-S tunneling) or 



both independent of whet 
X = +00 (S-N tunneling) [33|. Hence, W and 5 in the S-N tunneling case are the same as 
those in the case of N-S tunneling in Fig. [TOl 

The transmission probability for quasiparticle current Wj^ is shown in Fig. [TlTb). We 
flnd that Wj^ is enhanced to be greater than unity at low energies, due to the factor 
L(a2 + b^) = + {gnoY/E in Eq. ([771). 

When the Bogoliubov phonons propagate toward the S-N phase boundary. Fig. [T2r b) 
shows that atoms evaporate from the surface. We also flnd that Jq(x) changes near the 
surface and the supercurrent A Jg is induced, which flows toward the boundary. S{x) becomes 
positive around the surface of the superfluid region, reflecting the behavior of Jq(x) and 
/S.Js{x). This supercurrent AJs(x) is considered to be induced by the reflected Bogoliubov 
mode, which twists the condensate phase in the superfluid region. The fact of Rj + Wj > 1 
reflects that Jq{x) increases during the tunneling through the S-N phase boundary. 
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VI. CONCLUSIONS 



To summarize, we have investigated tunneling effects of Bogoliubov excitations at T = 0. 
We have extended our previous work to the case when the condensate densities are different 
on the left and right of the barrier. Within the frame work of the Bogoliubov theory, 
we have evaluated the transmission probability, phase shift as well as the energy flux and 
quasiparticle current carried by Bogoliubov excitations. We showed that, while the energy 
flux is conserved, the quasiparticle current is not conserved. The excess quasiparticle current 
is actually cancelled out by the counterflow of supercurrent, which is induced by the back- 
reaction effects of Bogohubov phonons on the condensate. In the case of a rectangular 
potential barrier, we directly showed that the induced supercurrent satisfles the Josephson 
relation with respect to the twisted phase by Bogoliubov phonons. When the condensate 
has different densities on the left and right of the barrier, the supercurrent is induced in 
the region far from the barrier potential. We also studied the tunneling of atoms from 
the normal region to the superfluid region, as well as the tunneling of excitations from the 
superfluid region to the normal region. In the former case, we showed that supercurrent is 
induced inside a condensate by injecting free atoms from outside. In the latter case, we found 
that atoms evaporate from the superfluid-normal state phase boundary, when Bogoliubov 
excitations propagate toward the surface of the superfluid region. We think these results can 
be of interest for the investigation of Bogoliubov mode and its connection to the superfluidity 
of BECs in ultracold atomic gases. 
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APPENDIX A: FORMALISM OF WEAKLY INTERACTING BOSE GASES 



m 



n this appendix, we summarize the formahsm of weakly interacting Bose gases developed 
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351]. We introduce approximations for inhomogeneous Bose condensates 



including the Bogoliubov approximation used in this paper. 

We consider an interacting Bose gas described by the Hamiltonian, 

V2 



K 



dr ip^ 



2^ + U{r) - fij ij{r) + ^ I dr ij\r)ij\r)tlj{r)ij{r), (Al) 



where ipii^) is the Bose field operator, fi is the chemical potential, and U{r) is an external 
potential. We assume a contact interaction between atoms g6{r — r') with the coupling 
constant g = Anas/m, where a,, > is the s-wave scattering length. 

In the Bose condensed phase, we divide the field operator '?/'(r) into the sum of the 
condensate wave function "^oi^r) = {ip{r)) and the fluctuation part 6ip as 



(A2) 



Substituting Eq. flA2p into Eq. flAip . we approximately evaluate the cubic and quartic terms 
with respect to Sip and Sip'' as 



(A3) 
(A4) 



where n(r) = {dip'^dijj) is the non-condensate density, and m{r) = {dipdip) is the so-called 
anomalous average 35|. In this mean-field approximation, Eq. (lAip reduces to 



K = K0 + K1 + K2, 



Ki 
K2 



9 



dr %T^o + ^ / dr |^o| , 



dr 
dr 



6ip^ (f + 2^(|^oP + ^)) + I ((^0 + m)6'^U^^ + h.c. 



(A5) 
(A6) 

(A7) 

(A8) 



Here, T = —^ + U (r) — yU. From the condition that the linear term in terms of Sip and Sip''' 
vanishes, we obtain the (generalized) Gross-Pitaevskii equation 2J, |26|, 



2m 



-7^ + f/(r)+^7(|^o|' + 2n) 



^0 + grh% = /i^c 



(A9) 
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The quadratic term K2 in Eq. flASp can be diagonalized by the Bogohubov transformation 



25| 



r) a'j-Vj\T)aj 



(AlO) 
(All) 



where cij is the creation operator of a Bogohubov excitation in the jth state, which obeys 
the bosonic commutation relations, 



ti 



0. 



(A12) 



Diagonahzation of -Ko is achieved when [uj{r)^Vj(r)) satisfy the following generalized Bo- 



gohubov equations 

2m 
2m 

Then, we have 
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35|: 



+ f/(r) +2(7(|^or + ^) -/i 
+ f/(r) + 2(7(|^o|' + ^)-/W 



V, - g{{%y+m*)u 



-EjVj. 



K = Kn 



dr \vj\^ 



(A13) 
(A14) 

(A15) 



The last term in Eq. (lAlSP is the so-called quantum depletion, describing the non-condensate 
due to the repulsive interaction between atoms. It remains finite even at T = 0, where 

We note that Eq. flA9p involves terms originating from excitations {2gh'^Q and gm'^'^. 
This reflects the fact that the condensate and excitations affect each other. In the Bogohubov 
approximation, both n and m are neglected, so that effects of Bogohubov excitations on 
the condensate are not taken into account. Equation flA9P without n and m is the ordinary 



static GP equation 
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26l |. while Eqs. (lAlSp and (]A14p without m and n are the Bogoliub 



ov 



coupled equations. The Bogohubov approximation is vahd at very low temperatures where 
n and m are very small. 

The approximation keeping both n and m is the Hartree-Fock-Bogoliubov approximation 



351]. This approximation is valid at finite temperatures where non-condensate fluctuation 
cannot be neglected. In this approximation, however, the excitation spectrum in a uniform 



system has an energy gap 



35| . This is inconsistent with the Hugenholtz-Pines theorem 36| , 
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which states that the excitation spectrum must be gapless in the BEC phase. Keeping n but 



neglecting m is referred to as the Popov approximation 



35j. This approximation is also 



considered to be valid at finite temperatures. Since it yields a gapless excitation spectrum, 
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it has been widely used in the study of BEC at finite temperatures js], 
properties of Bogoliubov excitations at finite temperatures have been also studied using the 
Popov approximation 



APPENDIX B: CONSERVATION LAWS FOR BOGOLIUBOV EXCITATIONS 

In this appendix, we discuss the conservation laws in terms of quasiparticle current and 
energy fiux associated with Bogoliubov excitations. In the Bogoliubov mean-field approxi- 
mation, the total number density n = (ij^'^ip) and total current density J = (l/m)Im(-?/'"''V^/') 
are, respectively, given by 

n = ris + ^ {riu^ + n^^) + ^n^,, (Bl) 

j j 

j 3 

In this appendix, the index for eigenstates j is explicitly written. Note that in our tunneling 
problem of Bogoliubov excitation, each eigenstate in Eqs. fIBip and flB2p corresponds to 
Bogoliubov excitation with energy E injected from a; = — cxDora; = cxD. Here, Us = |^oP 
describes the condensate density and 

J, = -Im(^*V^o) (B3) 

m 

is the supercurrent density carried by the condensate, n^^, tt,^,^, J„^, and J^. are, respectively, 
given by 



Uu, = 




(B4) 


n,. = 




(B5) 


Juj = 


— \miu*Vuj), 
m 


(B6) 


Jvj = 


— lm(v*Vvj). 
m •' 


(B7) 



The total number density n and total current density J satisfy the continuity equation 

dtn + V ■ J = 0. (B8) 
31 



Since the second terms in Eqs. (IBip and (lB2p describe the quasiparticle contributions, the 
quasiparticle density riqj and quasiparticle current Jqj are, respectively, given by 

riqj = riu^ +n^., (B9) 

Jq,j = Juj — Jvj- (BIO) 

Equations (1B9P and (IBlOp show that both the quasiparticle density nq .,■ and current Jq^- 
consist of two components originating from Uj and Vj. We note that the current density of 
f component appears as —Jvj in Eq- (1B10I) . Thus, in a uniform system, a Bogoliubov phonon 
is accompanied by two current components, Ju^ = {p/m)a? and —Jv^ = —{p/m)b'^, where 
a and b are given in Eq. (1301) . Hence, the f component flows in the oppose direction to the 

m. The 



u component. Indeed, these counterpropagating currents were recently observed 
last terms in Eqs. fIBip and flB2p describe effects of quantum depletion. 

To derive the continuity equation for quasiparticles, it is convenient to use the time- 



dependent Bogoliubov equations 



39| for {u{r,t),v{r,t)), 




*\2 



^rsdtl 1 = 1 :"^| I I, (Bll) 



where h = + f/(r) + 2g\'^o? - Equation (IBllD reduces to Eqs. and (HD in the 
stationary state, {u{r,t),v{r,t)) = e~^^^^{uj{r),Vj{r)). 

Using Eq. (IBlip . one obtains the continuity equations for Uj and Vj, as 

dtn^^+V-J^^ = ^, (B12) 
dtn,^-V-J,^ = ^, (B13) 

where 

Sj = -%Im {-^lupj) . (B14) 
Thus, the continuity equation for quasiparticles is given by 

dtuq, + V ■ Jq,, = S, . (B15) 



In Eq. ( IBlSp . Sj works as a source term. This means that the total number of quasiparticles 



is not conserved when Sj 7^ 0. In a uniform system, one finds 5*^ = 0, so that the number of 
quasiparticles is conserved. On the other hand, since the source term Sj is finite near the 
potential barrier in our tunneling problem, the number of quasiparticles is not conserved. 
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We next consider the energy flux. For this purpose, we deflne the energy density operator 
p as 

p = ^p^fip + (B16) 

where T is deflned below Eq. (]A8|) . Using the Heisenberg equation idti^ = T-i^ + g-ip'^ 'ijj'ip ^ one 
obtains the continuity equation for energy density p as, 



dtp + V ■Q = Q. 



(B17) 



Here, Q is the energy flux operator, deflned by 



Q = — {y^^)m + g^^^^] -h.c. 
2m L V / . 



-—Re 

m 



(B18) 



Substituting Eq. ( ]A2p into p = (p) and retaining terms up to 0(5-?/'^), we obtain 



9 



(B19) 



In obtaining Eq. (lB19p . we have used {6ij) = 0. Using Eqs. (lAlOl) . (1A13I) . and (1A14I) . we 
obtain 



j 



where 



Po = ^o^^o + ||^o|' 



(B20) 



(B21) 



is the condensate energy density. Since the energy density p is a real quantity, the last term 
in Eq. flB20p must vanish, which gives 



(B22) 



The energy flux Q = (Q) can be also calculated in the same manner. The result is 

3 j 

Here, 



2m - 



f^o + ^7(|^or + 2n)^o + gm%) (V^S) - c.c. 



(B24) 
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is interpreted as the energy flux carried by the condensate. Actually, Qo identically vanishes 
when \E'o satisfles the (generalized) GP equation. 

The second terms in Eqs. (lB20p and (1B23P describe the quasiparticle contributions. Thus, 
the energy density for quasiparticles Pqj and energy flux for quasiparticles Q^j are, respec- 
tively, given by 

Pqj = Ej{nu^ - Uy^), (B25) 
Qq,, =E,(J„^+J,J. (B26) 
Equation (IB25P shows that the v component has a negative energy density —Ejn^.. In 



Eq. (lB26p . the v component appears as +EjJ^., which is in contrast to J^j in Eq. (IBlOp . 
where the v component appears as —Jvj- This is because the v component has a negative 
energy —Ej and counterpropagating current density —Jvj ■ In contrast to the nonconserved 
quasiparticle number density in Eq. (IBlSp . the continuity equation with respect to the energy 
density has no source term, as 

dtPcj + V ■ Qqj = 0. (B27) 

Namely, Qq.j is conserved everywhere in the stationary state. 

To examine the origin of the source term Sj in Eq. (IBlSP , it is convenient to consider the 
divergence of Eq. flB2p in the stationary state, 

V ■ J = V ■ J, + ^Sj {a]aj). (B28) 
i 

In obtaining Eq. flB28D . we have used Eqs. flBT2|) - flBT4D and flB22D . The static GP Eq. ^ 
guarantees the conservation of the supercurrent (V ■ Js = 0), so that Eq. (1B28P contradicts 
with the conservation of the total current J obtained from Eq. ( IBSp . unless the last term in 
Eq. flB28P vanishes identically. 

This inconsistency arises because effects of quasiparticles on condensates (back-reaction 
effect) are completely neglected in the Bogoliubov approximation. This problem can be 
solved by including quasiparticle contribution to the condensate in the GP equation as 

("S + U{r) + g\moA ^0 - 2gY,UjV* = p^o • (B29) 

In this modified GP equation, the last term on the left-hand side originates from the anoma- 
lous average m in Eq. flA9p as 

g{5i,5i,)% = -2g J2 Ujv*{a]a,)% - gj^ ^j^jK (B30) 
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Using Eq. (lB29p . the conservation of the supercurrent (V ■ Js = 0) is modified to be 

V- Js = -^S',(d]«,). (B31) 

Substituting Eq. flBSip into Eq. flB28p . we obtain the expected conservation of the total 
current V ■ J = 0. In the case of the one-dimensional model we are using in this paper, 
when we integrate Eq. flB3ip in terms of x from — oo to x, we obtain Eq. 
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